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Electron-spin nitrogen-vacancy color centers in diamond are a natural candidate to act as a 
quantum memory for superconducting qubits because of their large collective coupling and long 
coherence times. We report here the first demonstration of strong coupling and coherent exchange 
of a single quantum of energy between a flux-qubit and an ensemble of nitrogen-vacancy color 
centers. 



During the last decade, research into superconduct- 
ing quantum bits (qubits) based on Josephson junctions 
has made rapid progress [l]. Many foundational exper- 
iments have been performed [2HTQ] and superconduct- 
ing qubits are now considered one of the most promising 
systems for quantum information processing. However, 
the experimentally reported coherence times are likely to 
be insufficient for future large scale quantum computa- 
tion. A natural solution is a dedicated engineered quan- 
tum memory based on atomic and molecular systems. 
Since macroscopic quantum coherence was first demon- 
strated in Josephson junction circuits [2 , the question 
of whether or not coherent quantum coupling between 
a single macroscopic artificial atom and an ensemble of 
natural atoms or molecules is possible has attracted sig- 
nificant attention pTtfM] . In this letter we present for the 
first time evidence of coherent strong coupling between a 
single macroscopic superconducting artificial atom (a flux 
qubit) and an ensemble of electron-spin nitrogen- vacancy 
color centers (NV~ centers) in diamond. Furthermore, 
we have observed coherent exchange of a single quan- 
tum of energy between a flux qubit and a macroscopic 
ensemble consisting of - 3 X 10^ of NV- centers. This 
provides a foundation for future quantum memories and 
hybrid devices coupling microwave and optical systems. 

With the early successes of single atom quantum state 
manipulation [1^1, research in quantum information pro- 
cessing with atomic and solid-state systems has pro- 
gressed largely in a separate fashion. In recent years, 
significant effort has been devoted to coupling atomic 
and molecular system to solid-state qubits to form hy- 
brid quantum devices pTHTSj . Hybrid devices involving 
the integration of an atomic system with a superconduct- 
ing transmission line resonator have been realized [16- 
[T8] . Such schemes have the potential to couple super- 
conducting solid-state qubits to optical fields via atomic 
systems, thus enabling quantum media conversion. The 



coupling strength g of an individual atomic system to 
one electromagnetic mode in a resonator circuit is usually 
too small for the coherent exchange of quantum informa- 
tion. However, the coupling strength of an ensemble of N 
such atomic systems will be enhanced by a factor of y/N 
[19], allowing one to reach the strong coupling regime 
{gyfN ^ A>:,7, where ti and 7 are the damping rates of 
resonator circuit and atomic system). 

Of the many possible hybrid systems, coupling a flux- 
qubit to an (NV~) center in diamond is particularly ap- 
pealing. Firstly, the magnetic coupling between a flux- 
qubit and a single NV~ center can be three orders of mag- 
nitude larger than that associated with a superconduct- 
ing transmission line resonator [14l. Second, the ground 
state of an NV~ center is a spin 1 triplet due to its C^y 
symmetry (Figure lb). The S = 1 spin triplet \ms = 0) 
state is separated by 2.88 GHz from the near degenerate 
excited states \ms = il) under zero magnetic field (Fig- 
ure Ic). This energy separation is ideal for a gap tunable 
flux-qubit to be brought on and off resonance with it. 

In this Letter, we report on the first observation of 
vacuum Rabi oscillations between a flux-qubit and an 
ensemble of approximately three million NV~ centers in 
diamond. This demonstrates strong coherent coupling 
between two dissimilar quantum systems with an effec- 
tive collective coupling constant of ^fens ^ 70 MHz. 

We begin by describing our experimental setup as de- 
picted in (Figure 1). An NV~ diamond sample was 
prepared by ion implantation of 12 (^2+ jqq j^^y 
der high vacuum into high-pressure, high-temperature 
(HPHT)-synthesized type lb (001) surface orientation 
single crystal diamond. The 12^2+ with a dose 

condition of 3 x 10^^ cm~^, were stopped at a depth of 
600^^00 "^^^^ generated on the order of 5 x 10^^ cm~^ 
vacancies over a depth of ^ 0.7/im. After implantation, 
the crystals were annealed at 900° C under vacuum for 3 
hours. This high dose carbon implantation method en- 
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FIG. 1. Experimental set-up of an NV-diamond sam- 
ple attached to a flux-qubit system. (a) A diamond 
crystal is glued on top of a flux-qubit and its superconduct- 
ing circuits (under the red box) with the diamonds ^^C im- 
planted (001) surface facing the chip. The distance between 
the flux-qubit and surface of the diamond crystal is carefully 
adjusted to be less than a micrometer using 100 nm height 
mesa structure on the diamond surface, a circle adjacent to 
the red square, and the optical interference pattern (Newton's 
ring), (b) A sketch of a NV color center of diamond with its 
vacancy (V) and nitrogen (N) atom, as well neighboring car- 
bon atoms. Four equivalent NV-axis exist depicted in purple 
color dangling bonds, all making the same angle with [001] 
direction to which the magnetic field generated by the flux- 
qubit points, (c) Energy diagram of the NV center, with the 
spin triplet \ms — 0) ground state separated by 2.88 GHz 
from the degenerated \ms = ±1) excited states under zero 
magnetic field 28 . (d) Optical micrograph and the circuit 
scheme of the aluminum made flux-qubit, the magnified view 
of the chip under the red box region shown in Figla. The 
central M shaped circuit contains a flux-qubit and a SQUID 
detector. Two high-bandwidth (20 GHz) MW-control lines 
located both sides of the qubit circuit, (e) The H-shaped 
gap tunable flux-qubit and the edge-shared SQUID used as a 
switching qubit state detector. The flux-qubit contains two 
loops, the main loop (blue) and the a-control loop (magenta) 
which controls tunneling energy gap of the flux-qubit. Magni- 
fied view of Josephson junctions are also shown. The mutual 
inductance of control line-1 to the a- loop and main loop are 
90 fH and 256 fH, and those of control line- 2 are 0.5 fH and 
549 fH. The magnetic flux penetrating these two loops can be 
controlled in ns time scale by applying synchronized current 
pulses to these control lines in situ. 



hances the yield of generated NV~ centers [20^ . Photolu- 
minescence (PL) optical spectroscopy (shown in Figure 
2) established that NV~ centers were generated with a 
density of ~ 1.1 x 10^^ cm~^ over a 1/im depth. We can 
describe the ground state of a single NV~ center by the 
Hamiltonian [21 , 

i^Nv = hDSl + hE{Sl - SD + /i^nvMbB • S, (1) 

where Sx^Sy^Sz are the usual Pauli spin 1 operators, D 



the zero-field splitting (2.878 GHz), E the strain-induced 
splitting (<1 MHz), the N-V Lande factor g^v = 2, and 
= MHz/mT. The last term represents the Zeemen 
splitting, which is negligible in our case as the magnetic 
field applied perpendicular to the surface of the chip to 
prepare the flux qubit is less than 0.1 mT. 
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FIG. 2. Photoluminescence spectra of (a) ensembles of 
color centers in the highly carbon implanted sample and (b) 
a single NV- center in pure diamond at room temperature. 
Their signal intensities are normalized for the comparison of 
the spectra. In the spectrum of the single NV center, the con- 
tributions of phonon Raman scattering from bulk diamond at 
573 and between 600 and 620 nm are subtracted ^29^ . The zero 
phonon line of NV- at 637 nm[291 130j is clearly observed in 
both spectra. Broad spectrum of phonon replicas is also very 
similar to each other and to the reported ones 29, 30 .These 
indicate that the NV- center is produced as a major color 
center in the highly carbon implanted sample. The signal in- 
tensity of the ensemble is about 6.5 x 10^ times stronger than 
that of the single one. From this result, the concentration of 
the NV- center in the highly carbon implanted sample was 
estimated to be 1.1 x 10^^ cm~^. 

A diamond crystal was glued on top of the supercon- 
ducting circuit with the ^^C implanted surface facing the 
fiux-qubit (Figure la). We used a gap tunable flux-qubit 
[22] [22 (Figure Id, le) where the smallest junction of 
the three Josephson junction qubit is replaced by a low 
inductance superconducting quantum interference device 
(dc-SQUID) loop (the magenta loop in Figure le). The 
flux qubit - NV ensemble coupled system is measured by 
the qubit state using a built-in dc-SQUID (the biggest 
loop in Fig.le sharing edges with the flux-qubit) which is 
inductively coupled to the qubit. When biasing the main 
loop close to half a fiux quantum, the device is an effec- 
tive two-level system [24] described by the Hamiltonian 

i^qb = ^(Aa^ + e(j^), (2) 

which is given in the basis of clockwise and counter- 
clockwise currents. Here, the Pauli spin | matri- 
ces, he = 2/p(<l>ex — 3^o/2) is the energy bias (/p ~ 300 



3 



nA is the persistent current in the qubit, is the ex- 
ternal flux threading the qubit loop, and <l>o = h/2e 
is the flux quantum), and A is the tunnel splitting. 
The energy splitting of the gap tunable flux-qubit is 
hF = h\/ + where e and A can be controlled in- 
dependently by the external magnetic flux threading the 
two loops. This type of flux-qubit can be tuned into res- 
onance with an NV~ ensemble in-situ at a base temper- 
ature of ~12 mK while keeping the qubit at its optimum 
flux bias (degeneracy point). The total Hamiltonian of 
the coupled system is 

U= ^(Acr^ + ecr^) 

i 

+ \Y.3i'yzS^,u (3) 

where i runs over the NV~ centers which couple to the 
flux qubit. The corresponding coupling constant can be 
estimated using the Biot-Savart law at gi ~ 8.8 kHz. In 
our situation here, the | states of the NV~ electronic 
spin are near degenerate and so our flux-qubit couples to 
both the |0). ^ and |0). ^ transitions. This 

results in an effective coupling constant \f2gi larger that 
generally anticipated. 

From the spectroscopic measurements, a clear anti- 
crossing was observed (Figure 3a) near the degeneracy 
point of the flux-qubit, while no gap was observed in 
the same flux-qubit prior to the mounting of the ensem- 
ble (inset in Figure 3a). We also note a narrow reso- 
nance at 2.878 GHz of less than 1 MHz width near these 
anti-crossings. This can be ascribed to the near degener- 
ate excited states of the NV~ ensemble and so indicates 
a strain-induced zero-field splitting coefficient E of less 
than 1 MHz. From the fine scan spectroscopy shown in 
(Figure 3b), a vacuum Rabi splitting near ^ens ^ 70 MHz 
was clearly observed confirming strong coupling between 
the flux-qubit and the NV~ ensemble. Next from the 
measured vacuum Rabi splitting and our calculated value 
of Qi we can estimate the number of NV~ centers in the 
ensemble at A/" = gl^^/^g^ ^ 3.2 x 10^, where the factor of 
2 in the denominator is due to the two-fold degeneracy of 
the excited | states of an NV~ center. This estimate 
is consistent with the density of NV~ centers measured 
by PL spectroscopy in the whole sample (1.1 x 10^^ cm~^) 
multiplied by the volume of centers coupling to the flux 
qubit (area 40 jirn^x effective thickness 0.7 jj^m). The PL 
spectroscopy approach gives the number of coupled cen- 
ters as ^ 3.1 X 10'^. 

Next, we investigated the dynamics of our system in 
the time domain using a similar measurement cycle to 
that performed in qubit-LC resonator coupled systems 
[25] . We first excited the flux-qubit and then brought 
it into resonance with the NV~ ensemble. Single energy 
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FIG. 3. Energy spectrum of the flux-qubit coupled 
to ensemble of NV-centers. a) Resonant frequencies indi- 
cated by peaks in the SQUID detector switching probability 
when a 500 ns long microwave pulse excites the system be- 
fore the readout pulse. Data are represented as a function 
of the external magnetic flux through the effective qubit area 
i^ex = + where $m and are the flux through 

the qubit main-loop and a-loop, respectively). Inset: reso- 
nant frequency spectrum of the same flux-qubit over the same 
region without mounting the NV-diamond crystal, (b) Mag- 
nified view of the white dotted box region in (a). A vac- 
uum Rabi splitting as large as 70 MHz is clearly observed. 
Since the qubit phase relaxation time measured by spin-echo 
jiecho ^ Q 25 jjg^ ii^Q strong coupling condition is satisfied by 
the a/TV enhancement in ^. A/" 3 x 10^ is the estimated 
number of electron spins in the ensemble strongly coupled to 
the flux-qubit. 



quantum exchange between the flux-qubit and NV~ en- 
semble at resonance manifests itself as the vacuum Rabi 
oscillations 

|l)qb|0)ens-|0)qb|l)ens (4) 

where |l)ens = ^ Ei ^+,i|00 • • • 0), with ^+,, = 
|l)i(0|i + I — l)i(0|i being the raising operator of the i-th 
NV~ spin to both the | zb 1) states. (Figure 4a) clearly 
shows vacuum Rabi oscillations between the flux-qubit 
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FIG. 4. Vacuum Rabi oscillation of the flux-qubit 
NV-ensemble coupled system, (a) Damped oscillation 
(blue curve) of an initially excited flux coupled to an NV- 
ensemble with the results from a phenomenological model 
shown as a thin red curve. Inset, schematic measurement se- 
quence for the vacuum Rabi measurement. An Initial ground 
state |0)qubit|0)ens is prepared by letting the system to re- 
lax longer than 800 /is under the base temperature (T=12 
mK) of the dilution refrigerator at the optimal readout flux 
bias. Immediately following this a vr-pulse is applied to the 
flux qubit resulting in the state | l)qubit |0)ens- The system 
is then brought into resonance non-adiabatically by a flux- 
bias shift current pulse through the control line 2. On reso- 
nance vacuum Rabi oscillation takes place for a given time; 
|l)qubit|0)ens ^ 1 0) qubit | l)ens • Finally the SQUID detector 
reads the qubit state as a function of the time keeping the 
system on resonance. We repeat these measurement typically 
2 X 10^ times under the same condition to obtain good statis- 
tics, (b) 2D plot of the SQUID detector switching probability 
as a function of both the detuning by flux bias shift and the 
time keeping the flux-qubit NV-ensemble system at a given 
detuning. The white broken line corresponds to the switching 
probability shown in a). 



and ensemble of electronic spins at the 2.878 GHz res- 
onance. The decay time of the oscillations however is 
approximately 20 ns. This is much shorter than the re- 
laxation time of both the fiux-qubit (Ti^qb ~ 150 ns) 



and the NV~ ensemble (Ti^nv ^ 10//s). As we tune the 
flux-qubit away from the 2.878 GHz resonance, the decay 
time associated with vacuum Rabi measurement becomes 
significantly longer (Figure 4b). From these results, one 
must conclude that a source of strong dephasing of un- 
known origin exists in the system near resonance. There 
are several likely sources. The most probable is a large 
electron spin ^ bath from the PI (nitrogen atom sub- 
stituting a carbon atom) centers present in our HPHT 
Ib-type diamond crystal used to prepare the NV~ en- 
semble. In our situation, where there is no external mag- 
netic field, the NV~ centers and PI centers naturally 
couple [14 . Hanson et. al [26 have shown an enhanced 
decay may result. The PI center issue can be eliminated 
to a large extent by applying an external magnetic field 
to split the |±1) NV" states. A 1 mT field could split 
these by approximately 15 MHz, detuning the PI centers 
and thus significantly improving the dephasing time of 
the coupled system. We can also decrease the number of 
PI centers in the sample (from 100 ppm to 1 ppm) by us- 
ing different synthesized diamond crystals. In addition, 
by using non HPHT Ib-type crystals we can remove the 
effect of other natural defects that may be present. Fi- 
nally there is also a strong hyperfine interaction ('^100 
MHz) between the NV~ electron spin and ^^C nuclear 
spins. Without the nuclear spins being initially polar- 
ized, unwanted dephasing will result. By polarizing the 
nuclear spins this source of dephasing can be removed. 
This should allow us to observe vacuum Rabi oscillations 
where we are limited by T2 of the flux-qubit. 

In conclusion, we have experimentally demonstrated 
strong coherent coupling between a flux-qubit and an en- 
semble of nitrogen- vacancy color centers in single crystal 
diamond. Furthermore, we have observed, via vacuum 
Rabi oscillations, the coherent exchange (transfer) of a 
single quantum of energy. This is the first step towards 
the realization of a long lived quantum memory for con- 
densed matter systems with an additional potential fu- 
ture application as an interface between the microwave 
and optical domains. 
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